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Abstract. Vector algebra is a powerful and needful tool for Physics but 
unfortunately, due to lack of mathematical skills, it becomes misleading for first 
undergraduate courses of science and engineering studies. Standard vector identities 
are usually proved using Cartesian components or geometrical arguments, accordingly. 
Instead, this work presents a new teaching strategy in order to derive symbolically 
vector identities without analytical expansions in components, cither explicitly or using 
indicial notation. This strategy is mainly based on the correspondence between three- 
dimensional vectors and skew-symmetric second-rank tensors. Hence, the derivations 
are performed from skew tensors and dyadic products, rather than cross products. 
Some examples of skew-symmetric tensors in Physics are illustrated. 
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1. Introduction 

Vector analysis [T] plays a key role in many branches of Physics: Mechanics, Fluid 
dynamics, Electromagnetism theory.., because it is a powerful mathematical tool that 
can express physical laws in invariant forms. Hence, learning of vector skills must be a 
priority goal for science and engineering students of undergraduate courses [2]. However, 
understanding of vectors often becomes intricate [3] due to the underlying mathematics, 
which can even hide the meaning of the involved physical quantities [1] . Common pitfalls 
are originated by the lack of mathematical resources for deriving vector identities. 

In undergraduate physics courses, the standard identities of introductory vector 
algebra are mostly proved either from geometrical arguments [5] or analytically using 
rectangular Cartesian components, and at best from the indicial notation PQE]- Unlike 
the analytical proofs, geometrical derivations are performed regardless of the coordinate 
system. Instead, the demonstrations based on the indicial notation are more elegant 
and compact although they require to handle complex symbolic expressions, without 
any physical insight into the problem at hand. 

We present an alternate approach of vector identity derivation based on the use of 
tensors and dyadic products rather than cross products. Tensor algebra using matrix 
format [7] become less cumbersome than indicial notation and further, the operations 
involving second-order tensors are readily understood as transformations of vectors. 

Hereafter, only for illustrative purposes, just first- and second-rank Cartesian 
tensors are considered, i.e. the three-dimensional space is Euclidean. Hence, the 
contravariant and covariant components are identical to one another because the metric 
tensor and conjugate metric tensor are equal to the identity matrix. Nevertheless, 
the derivations compiled in this text are equally valid for other metrices with minor 
modifications. 

2. Dyadics 

Aside from the well-known dot product (particular case of the inner product), a dyadic 
is formed by the outer or direct product of two vectors. The dyadic between the vectors 
a and b produces the following second-order tensor [7] of nine components: 



with i, j = 1, 2, 3 and where and bj are the respective Cartesian components of both 
operating vectors. Unlike the inner product or contraction, symbolized by a point, and 
the double inner product, symbolized by colon, no specific symbol is employed for the 
dyadic product. 

Since an arbitrary vector can be expressed as a linear combination of the unit vector 
basis {c-i}-^^ 2 3' arbitrary dyadic can be written into components from the concerning 
unit dyads {eiCj} . .^-^ 2 3 follows: 






(2) 
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where the summation convention is in effect for the repeated indices [T]. If the unit 
vectors Cj are mutually orthogonal, a special dyadic called the identical dyadic arises: 

1 = CiCj (3) 

where the summation convention is again invoked. This quantity is the second-order 
identity tensor of three-dimensional space. 

The inner product can be applied between vectors and second-order tensors as well, 
like a matrix product keeping their own properties. Thus, dyadics hold the following 
properties (derivation not shown): 

• ab = (j)3^ 

• (m) ■ 6 = (^a -bj c 

ab] = (c - a) b 



abj ■ (^c(tj = ' (^d(tj 



where the superscript t stands for the matrix transpose. Note that even though the 
vector transpose is represented by a 1 x 3 matrix instead of the conventional 3x1 matrix, 
the vector after transposition remains identical, i.e. a = [of . By default, vectors at 
left-hand side in an inner product are transpose. 

Although it is not used in this paper, the trace of ab, i.e. the sum of their diagonal 
components, is indeed the concerning dot product: 



trace yabj = a ■ b 

In fact, the trace of ab can be expressed in terms of double inner product as |a6 : 1. 
3. Skew-symmetric tensor associated to a vector 

In vector algebra [8], the skew-symmetric tensor of rank two associated to a vector 
a is defined by: 

{^a)ij = -Sijkau (4) 

where eijk stands for the Levi-Civita symbol also referred to as e-permutation 
symbol, and where all indices have the range 1, 2, 3. The index k is the dummy 
summation index according to the summation convention. The epsilon symbol Sijk 
holds the following rules: 

• £^123 = ^231 = ^321 = 1 

• ^123 = —£^213 = —£^132 

• ^ijk = 0, otherwise 
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There is an additional relation known as epsilon-delta identity: 

where 6ij is the Kronecker delta (^j-component of the second-order identity tensor) 
and the summation is performed over the i index. Indeed, the epsilon symbol and 
the Kronecker delta are both numerical tensors which have fixed components in every 
coordinate system. As the identity tensor,!, can be generated from the summation of 
the unit dyads ([3]) built by any orthonormal vector basis {ej}j^^2 3; epsilon symbol 
can be accordingly found from the following triple scalar product: 

where the cross product is symbolized by x. From the anti-cyclic rule of Sijk and the 
definition (jll), it is straightforwardly shown that the tensor Qg is anti- symmetric: 

oi = -n, (6) 

and this can be readily illustrated from the matrix form of fig: 
^ —as 02 ^ 

^a= «3 -ai 

y -a2 ai J 

Also, a is called the (Hodge) dual vector of the skew-symmetric tensor Og. Hence, for 
instance, the magnetic field tensor in Electrodynamics [U] is indeed the skew-symmetric 
tensor associated to the magnetic field vector. 

The Levi-Civita symbol also appears in the definition of the cross product of a and 

6PP: 



ax b) = Eijkajbk (7) 

/ i 

then, from the definition (jll) and the anti-cyclic rule of the epsilon symbol, it is possible 
rewritten the cross product in terms of the concerning skew-symmetric tensor (jl]) as: 



axbj^ = {n,)^,bk (8) 
or in vector notation as: 

a X 6 = fig ■ 6 (9) 

A cross product typically returns a (true) vector or polar vector. More exactly, the 
cross product ([H]) is a vector if either a or 6 (but not both) are pseudovectors. Otherwise, 
a X 6 is a pseudovector |10j . Then, it is worthy to mention that the tensor will be 
a relative tensor or pseudotensor if the vector a is axial and otherwise, it will be an 
absolute tensor if the vector a is polar. 

In addition to the skew-symmetry ([U]), the tensor fig holds the following properties 
(derivation not shown): 
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• Cls -3=0 

• ilg ■ a = b ■ fig: 

• fls • rig = 6a — (^a -bjl 

• ^dxb = ba - ab = fig ■ flj; - fl^ ■ Q,^ 

where a is a scalar. These properties can be straightforwardly proved using index 
notation and the above-mentioned rules of the Levi-Civita symbol. Thus, the epsilon- 
delta identity ([5]) draws to the last two properties, which are very helpful for the 
derivations compiled in section |H In particular, these other properties are also very 
useful: 

• Qs ■ b = -b ■ fig 

• nl = ee-l 

• nl = - fie 

where e is a vector of unit length. Due to Eq. ([3]), f2| is equal to the second-order 
identity tensor of the two-dimensional space (plane) with normal unit e. 

4. Standard vector identities 

Next, the most useful vector identities are demostrated from the concerning dyadics 
([T]) and skew-symmetric tensors (jlj). The above- listed properties, the associative rule of 
matrix product and the matrix transposition rules are used accordingly. 

• Cyclic permutation of the scalar triple product: 

(^a X bj ■ c= {a - fig) ■ c = a ■ (fig ■ c) = a ■ ^6 x 

= {b- fl^n) ■ c = b- {fl^a ■c) = b-{cxa) (10) 

From these identities, the orthogonality between a x b and each vector can be 
readily illustrated: 



dxbj-a=0 (11) 
Vector triple product expansion (or Lagrange's formula): 
ax \ b X cj = a ■ fig^? = S ■ ic^ — be 



= {a- c)b- (^a-bj c (12) 

also known as the acb minus abc rule. Using this identity, any vector can be 
expressed as linear combination of two mutually perpendicular vectors according 
to an arbitrary direction, e: 

a = (a ■ e) e + e X (axe) (13) 
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Furthermore, Eq. ( 1121) might be used directly in many below-mentioned identities, 
thereby it is one of the most important vector identities. 

Jacobi's identity: 

ax (^xcj+bx{cxa) + cx (^a x = a - ^2^^^ + b ■ ^sxa + 

+c-n.^g =a-(^S-bcj + 

+b ■ (dc — ca) + c- (ba — db) =0 (14) 



i.e. the sum of all the cyclic permutations of the vector double product comes to 
zero. 

Dot product of two cross products: 

(ax by (cx dj = (fig ■ ■ (cts-dj = b- (f^.g ■ O,-) ■ d 
= b ■ ((a ■ c) 1 — cd) ■ d 



b ■ (id ■ c) d— (d ■ ctj 

c)(b-dj - (d-dj (b-dj (15) 



[a ■ c 



This identity and the next one afford a simple means of deducing the formulae of 
Spherical Trigonometry. From Eq. (1151) . it is also derived the well-known identity: 

dxby = {d-d)(b-b^-(d-by (16) 

which geometrical interpretation is the Pythagorean theorem. 
Cross product of two cross products: 

d X b] X (^cx (tj = Og^g ■ (fl^ ■ ctj = (bd — dbj ■ (fts ■ ctj 

'bd^ ■ (Cts ■ (tj - (dbj ■ (fls ■ (tj 



'Qs- dj) b- (b- (n^-dj 
= (d- (cx (tj^b- (b- (cx dj^ d (17) 

(dxbjx (cxdj = (a ■ fig) • 0.^^-= (a ■ rig) ■ (dc- cdj 
= (a ■ fig) • (dcj - (a ■ Og) • (^cdj 
= (d- (n^-dj^ c - [d- {n^-c))d 
= (d- (bxdj^c- (d- (bx Syd (18) 



Other identities: 



d X b] X b] ■ d = { fi^xb ' bj ■ d = \^ ybd — db] ■ b] - a 
d ■b]b — \ b -b] d] ■ d 



-{dxbf (19) 
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The identity (1161) was invoked for this result. However, (^(^a x x ■ b = is 
readily found from Eq. (1111) . Likewise, in the next two expressions, the identity 
(JTTl) was again apphed accordingly: 

a X b^ X b^ X a = ri^'^^g^ xb' ^ ^ (^a x b^ — (^a x b^ bj ■ a 

a-b) (axb) (20) 



ax b] X b] X b = ^^Cc^xfe^xb ^ i b { a x b] — { a x b] b] ■ b 



-(b-b){axb) (21) 



5. Skew-symmetric tensors in Physics 

The substitution of physical pseudovectors (such as angular velocity or magnetic field) 
with skew-symmetric tensors (jlj) provides an alternate to cross product. This notation is 
much easier to work and allows to understand the vector operations in terms of rotations 
[TT] . In fact, an arbitrary vector a, which rotates an angle 9 about an axis along the 
unit vector e [l2], is expressed by: 

a* = (1 + ((1 - cos 6) Qg + sin ^1) ■ rjg) ■ a (22) 
5.1. Rotating systems 

A system rotates with constant angular velocity, uj, relative to a rest frame. The time 
variation of a unit vector e fixed to the rotating system [131 E| is given by: 

^ = uxe = na-e (23) 

where the tensor fl^ acts as a rotation operator. Also, if a point particle linked to 
the rotating system is moving at linear velocity, v, regarding to the rest frame, then it 
undergoes the following Coriolis' acceleration [131 [E]: 

acoriolis = 2UJ X V = n2Gj " V (24) 

where now the rotated vector is v. 



5.2. Rotation dynamics of rigid body motion 

A system of A^— point particles of mass rrii and position vector relative to a rest 
reference frame, describes a pure rotation with angular velocity u and acceleration a 

•Kinematics 

The linear velocity and acceleration of the i-particle are rewritten as: 

Vi = u X fi = il^ ■ fi (25) 

= dx ri + u X {u X n) = (O5 + fi^) ■ n (26) 
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•Inertia tensor 

The inertia tensor relative to the rest coordinate system is given by: 

N N 

I 4 ^ mj (rf 1 - fifi) = ^ m^rir- • (27) 

i=l 1=1 

•Steiner's theorem 

The inertia tensor relative to a second rest frame is: 



r = I - m}j fir-^ ■ n. 



(28) 



where fc is the position vector of the system center regarding to the initial rest coordinate 
system. 

•Inertia momentum with respect to an axis in the direction e 

N N 

h — ^^"^j (f^ X -Tj)^ = ^^"Zj (e ■ flfj ■ (r2_j?--e) = e ■ I ■ e (29) 

i=l i=l 

•Angular momentum 

N N N 

L = ^ruifi X {uj X fi) = ^rriiflf^- (r2_^vcJ) = (f2^-ri_^J ■ uj = I-u;(30) 

i=l i=l i=l 

•Rotation kinetic energy 

A 1 ^ 2 1 ^ 1 

Ek = - "^rriiiu} X r^f = - ^ (cl; ■ O^-J ■ (ri-r/t^) = -u-l-u (31) 

i=l i=l 

5.3. Electric quadrupole 

The quadrupole moment tensor of a system of point electric charges {qi}i^i ^ can be 
expressed as: 

Q 4 J] (3r>-;-r,2l) = 3 J] q^^l + 2 g.rf 1 (32) 

i=l i=l \i=l / 

This second-rank tensor is traceless. 

5.4. Vector field identities 

Vector field identities can be also derived using the skew-symmetric tensor associated 
to the differential vector operator V rather than the curl [15], as follows: 

Vxa = rJva (33) 

However, special care must be taken in tensor calculus because the order of elements is 
important: 

• (Va)* ^ aV 

• a ■ fiy 7^ — f^v ■ 
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• {aV)-b= (v-bja 

• (Va)*-6 = v) a 

Thus, the most relevant properties are: 

• riv ■ = (Va)* - (V ■ a) 1 

• ns-^v = (aV)*-l(a-V) 

• ■ (ns ■ 6) = (f^v • ^a) ■ b + {ns ■ n^f ■ b 

• V- = V) ■a + (fis- V)*-6 
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